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ABSTRACT 

A  Monte  Carlo  computer  simulation  was  developed  to  model 
hypothesized  electro-optic  laser  communication  systems 
operating  in  the  middle  ultraviolet  region  of  the  spectrum 
called  the  solar  blind.   By  assuming  various  source,  propa- 
gation, and  detector  characteristics  as  well  as  certain 
performance  parameters  it  is  possible  to  predict  the  effec- 
tive ranges  and  operating  characteristics  of  such  a  system. 
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I.   INTRODUCTION 

A.   HISTORY 

Historically  optical  communications  systems  were  devel- 
oped over  line-of-sight  paths  due  to  the  propagation  charac- 
teristics of  light.   In  the  United  States  Navy  line-of-sight 
communications  became  somewhat  neglected  with  the  advent  of 
uhf  radio  technology.   For  years  the  only  optical  communica- 
tion systems  in  the  navy  were  the  12  inch  (and  for  a  time 
24  inch)  hand  modulated  signal  lamps,  yardarm  blinkers,  flag 
hoist,  semaphore,  and  NANCY  (a  hand  modulated  infrared  sys- 
tem) .   Some  experimentation  was  conducted  in  the  late  1940's 
and  1950 's  using  cesium  lamps  operating  in  the  near  infrared 
(ir)  and  non-line-of-sight  or  over-the-horizon  (OTH)  commu- 
nication was  demonstrated.   Interest  waned,  however,  and 
progress  stood  still. 

Technology  has  not  stood  still  and  today  virtually  any 
radio  signal  can  be  monitored  and  the  transmitter  localized 
using  a  direction  finding  receiver.   While  communication 
systems  must  be  reliable  and  secure,  sometimes  it  is  impor- 
tant that  they  be  covert  (only  detectable  by  friendly 
forces).   Thus  armed  with  a  complete  set  of  radio  communica- 
tion equipment,  the  unit  commander  nevertheless  may  find 
his  forces  in  EMCON,  talking  to  each  other  via  flashing 
light. 
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B.  VISIBLE  LIGHT 

In  the  1960's  several  attempts  were  again  made  to  demon- 
strate the  feasibility  of  OTH  and  over-terrestrial-obstacles 
communication  using  atmospheric  scattering  [1-3] .   These 
efforts  utilized  visible  search  lights,  flash  lamps,  and 
lasers,  and  communication  at  night  was  successfully  demon- 
strated over  distances  up  to  fifty  kilometers. 

Visible  radiation  is  highly  attractive  because  of  the 
low  extinction  coefficients  (depending  on  weather  condi- 
tions) .   Unfortunately  it  is  seriously  degraded  during  day- 
light hours  due  to  the  background  radiation  from  the  sun 
which  tends  to  interfere  with  detection. 

C.  INFRARED  LIGHT 

Recently,  efforts  have  been  made  to  study  extended  line- 
of-sight  communications  utilizing  lasers  operating  in  the 
near  ir.   One  such  effort  utilized  a  1.06  y  laser  and 
employed  remotely  piloted  vehicles  (RPVs) .   Communications 
at  this  wavelength  appeared  quite  promising  and  ranges  from 
25  to  150  nautical  miles  have  been  predicted  [4] . 

In  this  region  of  the  spectrum,  however,  there  is 
noticeable  background  noise  due  to  scattering  of  incident 
solar  radiation  and  to  emission  by  atmospheric  particles 
heated  by  incident  radiation.   It  should  be  noted  that  the 
daylight  background  visible  radiance  is  due  primarily  to 
scattered  radiation  and  night  radiance  is  caused  mainly  by 
scattered  radiation  and  moon  light.   In  the  near  ir  daylight 
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background  radiance  is  due  mainly  to  atmospheric  emission, 
while  at  night  aurora  and  afterglow  are  the  predominate 
sources. 

D.   ULTRAVIOLET  LIGHT 

There  is  a  portion  of  the  spectrum,  however,  that  has 
been  studied  little  and  may  offer  much.   This  is  the  middle 
ultraviolet  (uv) .   Atmospheric  absorption  in  the  .2  to  .3  u 
region  is  caused  primarily  by  ozone  [5-6] .   The  relatively 
dense  layer  of  ozone  (at  approximately  22  km)  prevents  vir- 
tually all  of  the  radiation  in  the  220  to  280  nm  region  from 
reaching  the  surface  of  the  earth.   Thus  the  limiting  noise 
in  an  electro-optical  communication  system  operating  in  the 
middle  uv  would  be  self-noise  rather  than  atmospheric  noise. 
Another  critical  feature  to  be  considered  is  the  attenuation 
properties  of  the  intervening  atmosphere.   In  the  near  uv 
(wavelengths  between  .3  and  .4  y)  the  primary  attenuation 
mechanism  at  low  altitude  is  scattering.   This  can  be  divided 
into  Rayleigh  scattering  and  particulate  scattering  by 
atmospheric  aerosols.   Rayleigh  scattering,  inversely  pro- 
portional to  the  fourth  power  of  wavelength,  increases  quite 
rapidly  as  the  wavelengths  become  shorter.   In  practice,  it 
has  been  difficult  to  determine  the  particulate  scattering 
in  the  near  uv.   For  clear  sky  conditions,  however,  a  rea- 
sonable estimate  of  attenuation  can  be  made  by  assuming  that 
the  entire  attenuation  is  due  to  Rayleigh  scattering  [2] . 
In  the  middle  uv  particulate  scattering  is  very  sig- 
nificant.  At  high  altitudes  the  large  ozone  concentration 
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gives  rise  to  very  large  absorption  coefficients.   The  con- 
centration of  ozone  generally  decreases  as  the  altitude  is 
decreased.   Yet  even  at  sea  level  the  concentration  is  still 
high  enough  to  have  a  significant  effect.   These  propaga- 
tion characteristics  indicate  that  only  short  range  commun- 
ication systems  operating  in  the  middle  uv  are  likely  to  be 
useful.   Because  of  the  inherent  scattering  characteristics 
of  middle  uv  radiation,  along  with  the  minimal  background 
problem  in  the  solar  blind  region,  short  range  non-line-of- 
sight  communications  utilizing  broad  band  light  sources 
have  been  developed  [7-8] . 

Recently  a  number  of  reasonably  efficient  uv  lasers 
have  been  developed  with  energies  of  . 1  to  one  joule  per 
pulse  being  reported  [9-10] .   These  represent  a  new  genera- 
tion of  radiation  sources  to  be  considered  for  use  in  the 
solar  blind.   A  number  of  detectors  and  filters  have  been 
available  for  some  time.   It  was  these  considerations  that 
led  to  this  study. 

E.   PURPOSE 

The  purpose  of  this  effort  is  to  develop  an  analytical 
tool  for  estimating  the  ranges  for  which  optical  communica- 
tion is  possible  in  the  middle  uv. 
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II.   BACKGROUND 

A.   INTRODUCTION 

A  typical  communication  system  consists  of  a  transmit- 
ter, a  receiver,  and  the  intervening  space  between  them. 
The  characteristics  of  the  transmitter  and  the  receiver 
are  highly  dependent  on  each  other  and  on  the  properties 
of  the  propagating  medium.   Each  of  the  components  may, 
however,  be  treated  as  a  distinct  and  relatively  separate 
entity.   The  radiation  source,  or  transmitter,  must  be 
suitably  modulated  to  carry  information  and  must  have 
sufficient  power  to  overcome  the  attenuating  mechanisms  in 
the  medium  in  which  it  operates.   The  characteristics  of 
the  propagating  medium  are  not  well  defined  and  in  general 
must  be  approximated.   The  detector,  or  receiver,  must  have 
sufficient  sensitivity  to  detect  and  demodulate  the  atten- 
uated signal  that  is  embedded  in  background  noise.   Each  of 
these  components  is  present  in  an  electro-optical  communica- 
tion system  and  is  discussed  below.   It  will  be  seen  that 

1)  advances  in  sources  dictate  a  new  look  at  the  problem; 

2)  uncertainties  in  propagation  characteristics  of  the 
middle  uv  requires  careful  experimentation  to  gather  suitable 
information  such  that  a  proposed  system  may  be  modeled;  and 
that  3)   certain  aspects  of  detectors,  such  as  filters, 
critically  determine  the  performance  parameters  of  such  a 
system. 
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B.   SOURCES 

Ultraviolet  sources  suitable  for  communication  pur- 
poses include  omnidirectional  broad  band  lamps  such  as 
xenon  or  mercury  flash  lamps,  and  lasers.   Experiments 
utilizing  omnidirectional  sources  have  been  carried  out 
for  several  years  at  the  Aeronautical  Research  Associates 
of  Princeton  (ARAP) ,  Princeton,  New  Jersey  [7-8].   Uti- 
lizing a  pulsed  source  at  4  810  pps  and  an  average  power  level 
of  ten  to  thirty  watts ,  voice  communication  was  demonstrated 
over  a  distance  of  several  hundred  meters.   With  suitable 
modifications,  it  is  estimated  that  ranges  of  up  to  three 
kilometers  may  be  achieved  with  the  current  system.   With 
such  a  source,  the  conversion  efficiency  from  input  power 
to  output  power  is  about  0.1  percent. 

New  uv  laser  sources  reported  in  the  last  year  or  so 
appear  to  have  demonstrated  orders  of  magnitude  improvement 
in  efficiency.   Utilizing  rare  gas-halide  eximer  species, 
uv  output  with  working  efficiencies  of  up  to  ten  percent 
are  conceivable.   Practical  working  models  are  already 
commercially  available.   Wavelengths  in  the  range  of  200  to 
4  00  nm  have  been  reported,  indicating  that  various  wave- 
lengths in  the  middle  uv  may  be  attained  with  these  sources. 
The  power  output  of  these  devices  typically  runs  from  0.1 
joules  upwards,  depending  on  the  size  of  the  device.   The 
high  degree  of  collimation  of  the  output  suggests  that  the 
geometric  patterns  for  communication  will  be  highly  direc- 
tional.  This  is  one  aspect  of  the  problem  which  is  quite 
different  from  that  observed  for  omnidirectional  sources. 


15 


C.   PROPAGATION 

The  absorption  coefficient  for  middle  uv  radiation 
varies  by  orders  of  magnitude  as  a  function  of  wavelength. 
Scattering  by  particulates  at  these  wavelengths  is  more 
severe  than  in  the  visible  region  of  the  spectrum.   The 
combined  absorption  and  scattering  cross-sections  give  rise 
to  extinction  coefficients  which  have  been  measured  to  be 
in  the  range  of  0.1  to  five  per  kilometer  [11].   Further, 
variations  of  an  order  of  magnitude  in  extinction  coefficient 
over  a  spectral  range  of  less  than  fifty  nanometers  indicates 
extreme  sensitivity  to  wavelength,  due  largely  to  the 
absorption  bands  of  species  such  as  ozone  which  appear  in 
this  region. 

At  wavelengths  less  than  280  nm,  the  absorption  coef- 
ficient of  atmospheric  ozone  is  sufficiently  intense  as  to 
eliminate  the  large  majority  of  solar  radiation  from  the 
earth's  surface  [12].   This  is  largely  due  to  an  ozone 
concentration  centered  at  20  to  25  km  above  the  surface  of 
the  earth.   At  sea  level,  the  ozone  concentration  is  typi- 
cally an  order  of  magnitude  less  than  at  the  higher  alti- 
tudes.  This  gives  rise  to  reduced  solar  background  levels 
without  the  high  extinction  coefficients  found  at  higher 
altitudes.   Such  behavior  suggests  that  communication  in 
the  middle  uv  may  be  feasible  for  relatively  short  distances 
at  sea  level,  with  a  minimum  of  interference  from  solar 
background  radiation. 
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Further  major  considerations  for  middle  uv  propagation 
are  the  effects  due  to  multiple  scattering.   Although  the 
primary  light  beam  may  be  attenuated  by  an  extinction 
coefficient  which  is  readily  measured,  the  scattered  radia- 
tion, especially  after  multiple  scattering,  may  contribute 
considerable  intensity  to  the  radiation  density  at  an 
observer.   Shettle  and  Green  [12]  present  calculations  for 
solar  flux  as  a  function  of  the  sun's  angle  with  respect  to 
the  zenith.   Particularly  at  low  angles  of  incidence,  the 
diffuse  radiation  is  found  to  be  twenty  to  thirty  orders  of 
magnitude  more  intense  than  that  due  to  the  direct  solar 
radiation.   Such  figures  clearly  emphasize  the  importance 
of  properly  treating  the  effects  of  multiple  scattering 
when  taking  into  account  the  effects  of  particulates  in  the 
atmosphere. 

In  general  multiple  scattering  calculations  are  cumber- 
some to  handle  and  yield  limited  information.   Typically 
these  calculations  are  done  on  a  computer  using  either  an 
analytical  [13]  or  a  Monte  Carlo  [14-15]  appraoch.   The 
Monte  Carlo  calculations  are  relatively  easier  to  model, 
but  utilize  a  significant  amount  of  computer  time  and  only 
approximate  information  is  obtained.   Some  generally  useful 
results  have  been  presented  by  Bucher  [14-15]  again  indicating 
orders  of  magnitude  greater  radiation  densities  from  the 
multiple  scattering  effects  than  from  transmission  of  direct 
radiation.   The  available  results  from  such  calculations  do 
not  readily  answer  questions  such  as  1)   what  is  the 
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scattering  intensity  for  non-line-of-sight  communication 
given  an  obstacle  with  an  assumed  angle  between  transmitter- 
obstacle-receiver;  and  2)   what  is  the  actual  radiation 
density  as  a  function  of  field  of  view  and/or  position  of 
an  observer  located  off  from  the  axis  of  propagation  at  a 
given  distance  from  a  collimated  source  such  as  a  laser. 

D.   DETECTION/DEMODULATION 

The  information  available  from  a  pulsed  signal  depends 
on  the  pulse  width,  repetition  rate,  and  signal-to-noise 
ratio.   Considerations  such  as  these  are  included  in  papers 
by  Kennedy  [16-17]  in  which  communication  through  optical 

scattering  channels  is  discussed.   Considerations  leading 

i 
to  a  suitable  signal-to-noise  ratio  analysis  are  also  pre- 
sented by  Yarif  [18]  for  electro-optical  systems.   Whether 
the  information  is  pulsed  or  continuous  wave  greatly  affects 
the  reception  characteristics  of  the  signal.   The  background 
light  levels  coming  through  the  filters  utilized  before  the 
detector  critically  determine  the  signal-to-noise  ratio. 
The  characteristics  of  the  filter  and  detector  thus  are 
crucial  in  determining  the  information  available  in  a 
communication  system. 
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III.   STATEMENT  OF  THE  PROBLEM 

One  of  the  primary  constraints  on  the  development  and 
deployment  of  any  new  system  in  the  United  States  Navy  is 
the  limited  fiscal  resources  available  for  such  efforts. 
With  the  advent  of  new  technology  it  is  nevertheless  impos- 
sible to  utilize  the  technology  to  build  a  new  system 
without  first  having  a  mission  or  demonstrated  need  for  that 
system  and  second  having  a  well-defined  degree  of  confidence 
in  the  ability  of  that  system  to  meet  that  need. 

An  electro-optical  communication  system  utilizing  a 
laser  operating  in  the  middle  uv  is  such  a  system.   (It 
requires  little  or  no  imagination  to  visualize  the  advan- 
tages of  rapid  short  range,  possibly  non-line-of-sight , 
optical  communications  over  the  methods  presently  employed 
by  the  Fleet.)   While  it  is  beyond  the  scope  of  this 
project  to  develop  such  a  system,  it  is  within  the  scope 
of  this  effort  to  formulate  a  model  of  such  a  system  capable 
of  predicting  the  various  parameters  of  that  system. 

The  need  for  modeling  is  further  enhanced  by  the  limited 
amount  of  information  available  concerning  the  middle  uv. 
By  utilizing  available  information  and  by  characterizing 
source,  propagation,  filter  and  detector  characteristics,  it 
is  possible  to  formulate  a  model  that  will  yield  order  of 
magnitude  results  for  any  proposed  system.   In  the  process 
of  modeling  and  analysis  the  need  for  critical  experiments 
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which  would  allow  the  modeling  efforts  to  proceed  with 
a  greater  degree  of  accuracy  is  sure  to  be  revealed. 
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IV.   LASER  SOURCES 

The  past  two  or  three  years  have  seen  a  rapid  develop- 
ment of  lasers  operating  in  the  ultraviolet.   The  most 
notable  recent  advances  have  occurred  with  the  rare  gas- 
halide  lasers  in  which  an  excited  state  eximer  is  formed. 
With  species  such  as  XeF*,  the  ground  state  is  not  an  asso- 
ciated molecule,  thus  depopulation  of  the  lower  level  is 
not  a  barrier  to  population  inversion.   Further,  such  rare 
gas-halide  lasers  have  an  inherently  high  quantum  efficiency. 
Efficiencies  relative  to  the  energy  deposited  by  an  excitat- 
tion  electron  beam  are  as  high  as  ten  percent,  while  overall 
working  efficiencies  of  up  to  one  percent  have  been  reported 
[9,10,19].   With  higher  overall  efficiencies  anticipated, 
such  lasers  represent  an  attractive  source  of  radiation  for 
middle  uv  applications. 

The  observed  wavelengths  of  several  rare  gas-halide 
lasers  are  presented  in  Table  I:.:    It  should  be  noted  that 
the  laser  output  from  these  sources  may  be  Raman  shifted 
with  reasonable  high  conversion  efficiencies  [20] .   By 
passing  the  output  from  a  KrF  laser  through  ten  atmospheres 
of  hydrogen,  conversion  to  the  first  stokes  wavelength  of 
2  5  percent  and  to  the  second  stokes  wavelength  of  ten  percent 
was  observed  [20].   At  eighty  atmospheres  pressure,  more 
than  fifty  percent  conversion  efficiency  to  the  first 
stokes  wavelength  has  been  observed. 
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TABLE  I 

Characteristic  Wavelengths  of  Several  Rare 

Gas-Halide  Lasers  Along  with  Some  Stokes 
Shifted  Wavelengths  Observed  with  Hydrogen  Gas 


First  Stokes     Second  Stokes 
Eximer   Wavelength  (nm)   Wavelength  (nm)   Wavelength  (nm) 

279  318 


KrF 

249 

XeF 

351 

XeCl 

308 

XeBr 

282 

322 


The  practical  performance  of  rare  gas-halide  lasers  is 
illustrated  by  the  commercial  availability  of  a  model  which 
will  produce  0.1  joule  pulses  at  15  pps.   One  table  top 
model  is  reported  to  produce  1. 5  joule  pulses  with  a  125 
nanosecond  pulse  width.   The  development  of  these  lasers  is 
in  its  infancy.   Sources  with  higher  repetition  rates  and 
greater  reliability  are  to  be  expected. 

Due  to  the  rapidly  expanding  availability  of  middle  uv 
lasers  sources,  no  single  presently  available  source  is 
assumed,  but  rather  the  effort  in  the  succeeding  sections 
considers  the  possibility  of  utilizing  a  middle  uv  source 
at  representative  wavelengths  and  with  an  output  power 
consisting  of  milli joule  to  joule  pulses.   Thus  the  behavior 
of  any  unforseen  laser  or  coherent  beam  may  be  predicted. 
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V.   PROPAGATION  PROPERTIES  OF  THE  ATMOSPHERE 
AND  THE  MULTIPLE  SCATTERING  PROBLEM 

A.  INTRODUCTION 

A  laser  signal  propagating  through  the  atmosphere  is 
subjected  to  various  attenuation  mechanisms.   The  time 
spread  of  a  pulse  can  be  broadened  due  to  scattering,  while 
the  amplitude  can  be  decreased  due  to  absorption  and  scat- 
tering by  atmospheric  constituents.   The  inherent  direction- 
ality of  a  laser  source  indicates  that  the  position  of  the 
detector  relative  to  the  source  should  also  affect  the 
reception  of  the  signal.   Since  a  large  fraction  of  the  total 
extinction  process  is  due  to  scattering,  multiple  scattering 
should  be  an  important  consideration  affecting  the  propaga- 
tion characteristics  for  distances  greater  than  one  or  two 
extinction  lengths. 

A  complete  description  of  all  of  the  effects  of  the 
atmosphere  including  inhomogeneities  (clouds,  terrestrial 
objects,  surface  effects,  etc.)  is  beyond  the  scope  of  this 
paper.   A  relatively  concise  treatment  of  the  atmosphere 
and  the  multiple  scattering  problem  entails  characterizing 
the  atmosphere,  simulating  atmospheric  effects,  and  tabulating 
the  results  of  the  simulation. 

B.  CHARACTERISTICS  OF  THE  ATMOSPHERE 
1.   Attenuation 

In  order  to  characterize  the  propagation  of  optical 
information  in  the  atmosphere  it  was  first  necesary  to 
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determine  the  extinction,  absorption,  and  scattering  coef- 
ficients.  Table  II  was  constructed  based  primarily  on  the 
information  reported  by  Dunkelman  [21] .   The  Rayleigh 
scattering  coefficients  were  adapted  from  Penndorf  [22] . 

For  the  purpose  of  this  investigation  the  following 
definitions  and  relationships  apply: 
A    is  the  wavelength 
K    is  the  absorption  coefficient 
K    is  the  Rayleigh  scattering  coefficient 
Kqp   is  the  particulate  scattering  coefficient 
K    is  the  total  scattering  coefficient 
K    is  the  total  extinction  coefficient 

KS   =   KSR   +   KSP 

K    =   KS   +   KA 

R    is  the  ratio  of  K   to  K  or  (K  /K) 

R    is  the  ratio  of  K    to  Kg  or  (Ksp/Ks) 

From  Table  II  it  can  be  seen  that  K._  ,  K__,  and  K 

A   SR 

tend  to  decrease  as  wavelength  increases  in  the  middle  uv. 
The  relationship  of  K   and  consequently  K  ,  R  ,  and  R 
to  wavelength  is  not  so  easily  described. 

2.   Angular  Distribution  of  Scattered  Photons 

Complete  characterization  of  the  propagation  of 
optical  information  required  the  angular  distribution  or 
phase  function,  P(9),  of  the  scattered  photons  for  the 
wavelengths  of  interest.   Figure  1  [23]  was  the  primary 
source  for  phase  function  information. 
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ABSOLUTE     PHASE    FUNCTION 
vs    SCATTERING    ANGLE 


7/15/77 


x  220.0    NM 

A  248 . 5    NM 

O  265  .0    NM 

D  280.5     NM 


60  120 

Scattering  angle,  /3  (deg) 


FIGURE   1.      Absolute   Phase   Function   vs.    Scattering  Angle 
[23] 
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In  order  to  mathematically  model  this  information, 
the  Henyey-Greenstein  function,  Eq.  1,  and  some  of  its 
modifications  were  studied.  The  Neer-Sandri  function  [7] , 
Eq.  2,  was  a  logical  choice;  however,  the  function  used  by 
Zachor  [24],  Eq.  3,  was  chosen  for  simplicity  in  modifying 
the  author's  initial  phase  function  to  properly  account  for 
backscatter ing . 


P(6,g)   = 


_  i-g 


2  r 


4tt 


(l+g2-2g  cos  6)3/2 


(1) 


P(6,g)   = 


_  izai  r i 

471    .  (l+g2-2g 


cos  6) 


372 


g(3  cos  6-1) 


1        1     2     1        ■ 572 
2  cos  6rt  (1+g  -2g   cos  9  J  )  ' 


(2) 


where 


cos  eQ|   =  1/7 


P(6,g)   = 


_  i-g 


2  r 


4tt 


(l+g2-2g  cos  6)3/2 


f  0.5(3  cos  9-1) 


dV)3/2 


(3) 


Fitting  the  data  reported  by  Neer  [23]  with  Eq.  3 
resulted  in  Figures  2  and  3.   These  figures  represent  the 
absolute  and  normalized  phase  functions,  while  Table  III 
lists  the  coefficients,  utilized  for  this  simulation. 
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0  60  120 

SCATTERING  ANGLE,  9   (DEGREES) 


FIGURE  2.   Absolute  Phase  Function  vs.  Scattering  Angle 

Utilizing  a  Modified  Henyey-Greenstein  Function 
to  Approximate  the  Data  by  Neer,  et  al.  [23] 
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FIGURE  3.   Normalized  Phase  Function  vs.  Scattering  Angle 
Utilizing  a  Modified  Henyey-Greenstein  Function 
to  Approximate  the  Data  by  Neer ,  et  al .  [2  3] 


29 


H 
H 


9 


o 

•H 
-P 
rG 

rH 

I 

•H 
W 

0) 
X! 

-p 

c 

•H 

0) 

CO 
D 

co 
-P 
C 
(1) 
•H 
U 
•H 
M-l 

m 
cu 
o 
u 


p^ 


Pi 


co 


« 


03 

\ 

P* 
co 


CO 


« 


g 


^-» 

O^i 

rH 

VD 

1 

ID 

to     6 

O 

«       ^ 

• 

^-* 

r^ 

rH 

CO 

P^     1 

rH 

CO       g 

r- 

W         ^ 

• 

H 

CN 

Pi 

1 

CD 

CO 

S 

CO 

« 

44 

m 

44 


.-< 


in 


CO 


in 


CO 


m 
in 


as 


m 
o 

o 


o 

CN 

CN 

oo 


in 
H 


CTi 


CN 
t>- 


CN 
CO 


o 

m 


CN 

n 


in 
in 
o 

CN 


in 
oo 


CN 


in 

CN 

m 


oo 

CN 

o 


g 

o 

O 

O 

c 

m 

00 

O 

CN 

CN 

m 

cu 

13 

x; 

C 

p 

rfl 

cu 

in 

rH 

p- 

-H 

• 

X! 

o 

5 

cu 

u 

u 

•H 

<D 

p 

5 

CO 

•H 

-0 

g 

CU 

•H 

N 

-P 

•H 

a 

• 

rH 

0 

0 

•H 

•H 

-p 

cu 

-P 

3 

X! 

to 

•H 

M-l 

0 

g 

p 

•H 

T3 

to 

C 

u 

to 

na 

fa 

CU 

cu 

a 

0i 

a 

a 

p 

m 

rC 

<d 

0 

X! 

g 

S 

CO 

C 

CU 

CU 

g 

3 

o 

0 

iH 

m 

CO 

<d 

CN 

> 

cu 

u 

X} 

CU 

0 

X! 

IW 

0 

p 

CO 

p 

CO 

cu 

u 

0) 

3 

fO 

10 

rH 

cu 

rfl 

rd 

a 

U 

> 

rd 

CU 

CU 

cu 

A 

g 

M 

EH 

s 

X! 

P 

o 

• 

o 

rH 

>1 

m 

H 

■H 

rrj 

CU 

U 

> 

O 

u 

-H 

M-i 

0 

P 

Cu 

U 

CO 

CU 

CU 

a 

3 

CO 

rH 

CU 

<d 

U 

> 

30 


C.   DESCRIPTION  OF  THE  CHARACTERISTICS  OF  THE  SIMULATION 

The  purpose  of  this  simulation  was  to  provide  approxi- 
mate results  for  the  multiple  scattering  problem.   By 
assuming  various  source,  propagation,  and  detector  charac- 
teristics, a  user  could  in  principle  predict  the  effective- 
ness of  a  proposed  system. 

The  Monte  Carlo  method  described  in  Appendix  A  was  employed 
to  generate  information  detailing  characteristic  patterns  of 
photon  scattering  as  a  function  of  wavelength,  position  of 
an  observer,  pulse  spreading,  and  field  of  view  cf  the 
detector  being  employed.   The  observer  may  be  located  at  any 
arbitrary  angle  off  from  the  axis  of  propagation  and  at  any 
arbitrary  distance  from  the  source  (measured  in  extinction 
lengths) . 

Of  primary  interest  were  the  characteristic  patterns  of 
multiply  scattered  radiation.   The  simulation  provided  for 
obtaining  diagrams  of  photon  flux  at  a  given  position  as 
viewed  by  either  the  source  or  the  receiver.   In  addition, 
the  patterns  received  by  the  detector  were  capable  of  being 
limited  by  field  of  view  restrictions.   By  measuring  the 
photon  flux  at  a  given  position,  contours  of  equal  photon 
flux  were  generated.   These  contours  provided  a  simple 
means  of  visualizing  the  radiation  patterns. 

Detector  field  of  view  considerations  indicate  the 
field  of  view  required  to  receive  a  given  threshold  of 
photons  at  a  fixed  position  as  well  as  the  maximum  range  at 
a  given  angle  that  a  threshold  number  of  photons  can  be 
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detected  for  a  fixed  field  of  view.   Finally,  pulse  spreading 
information  due  to  the  time  delay  caused  by  multiple 
scattering  was  also  available. 

D.   CHARACTERISTIC  PATTERNS  OF  MULTIPLY  SCATTERED  RADIATION 

Figures  4  and  5  represent  the  pattern  of  photons 
received  by  a  180  degree  field  of  view  detector  located 
zero  to  nine  and  36  to  45  degrees  off  from  the  axis  of 
propagation  respectively  at  one  extinction  length  from  a 
300  nm  laser  source.   As  the  detector  approached  the  axis 
of  propagation  more  photons  were  received.   The  photons 
along  the  phi  prime  equal  to  zero  axis  represent  single 
scattered  radiation.   Many  of  the  photons  received  were, 
however,  multiply  scattered. 

Figures  4  and  5  are  similar  to  Figures  6  and  7  only 
in  the  latter  figures  the  receiver  is  located  five  extinc- 
tion lengths  from  the  source.   The  percentages  of  multiply 
scattered  photons  increases  greatly  as  the  receiver  is 
moved  either  off  from  the  axis  of  propagation  or  away  from 
the  source.   At  five  extinction  lengths,  virtually  none 
of  the  radiation  received  is  single  scattered.   This  is  a 
strong  indication  of  the  hazards  of  single  scattering  models 
in  the  middle  uv. 

Figures  8,  9,  and  10  represent  the  contours  of  equal 
photon  flux  for  250,  280,  and  300  nm  wavelengths.   The 
label  for  each  contour  is 
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180 


270 
RECEIVER   LOCATED   ONE   EXTINCTION   LENGTH 
FROM  THE   SOURCE.    9  =  0  TO  9   DEGREES 
A  =300  NM.  FIELD   OF   VIEW  =   180   DEGREES 


FIGURE  4.   Photon  Scattering  Diagram  (6»  vs  (J)')  as 
Viewed  by  a  Receiver  Located  0  Degrees 
from  the  Axis  of  Propagation. 
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270 
RECEIVER    LOCATED   ONE   EXTINCTION   LENGTH 
FROM   THE   SOURCE.    9  =  36   TO  45   DEGREES 
X  =  300  NM.   FIELD   OF   VIEW  =    180   DEGREES 


FIGURE  5.   Photon  Scattering  Diagram  (9'  vs  <J> ' )  as 
Viewed  by  a  Receiver  Located  9  Degrees 
from  the  Axis  of  Propagation. 
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270 
RECEIVER    LOCATED  FIVE   EXTINCTION   LENGTHS 
FROM  THE   SOURCE.    5=0T0   9   DEGREES 
X  =  300  NM.   FIELD   OF   VIEW  =  180   DEGREES 


FIGURE  6.  Photon  Scattering  Diagram  (9 
Viewed  by  a  Receiver  Located 
the  Axis  of  Propagation 


vs  <J> ' )  as 

)  Degrees  from 


35 


90 


180 


^60 

/               *    \           ^*00' 

4 

A^ 

+      ^^>s/ 

♦           \ 

+ 

4-     / 

\  0>^o\ 

1                   /  ^^^     *       S\ 
1                    1           ^Svs<    /        \ 

♦*"></ 

x/'^K         \ 

1                       \       *      ^ 

♦  ♦  i          1 

\                     ^s^                *  /^ 

\ 

+      yC.                / 

270 
RECEIVER    LOCATED   FIVE   EXTINCTION   LENGTHS 
FROM   THE   SOURCE.    0  =  36  TO  45  DEGREES 
X  =  300  NM.  FIELD   OF   VIEW  =  180  DEGREES 


FIGURE  7.   Photon  Scattering  Diagram  (9'  vs  (J)1)  as  Viewed 
by  a  Receiver  Located  9  Degrees  from  the  Axis 
of  Propagation 
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Relative  Flux  =   -  log. 


r         ■ 

1  N   2ttR   (cose, 
L  o  1 


rR   (COS0.,  -  cos0?) 


(4) 


2 

where  2ttR   (cos 6,  -  cos  02)  is  the  area  of  that  portion  of 

a  hemisphere  located  R  extinction  lengths  from  the  source 
subtended  by  0,  and  B-,  N  is  the  number  of  photons  pene- 
trating that  area,  and  N  is  the  number  of  photons  leaving 
the  source.   Each  successive  contour  of  photon  flux  represents 
a  decrease  of  one  order  of  magnitude  in  the  relative  photon 

flux.   Thus,  for  example,  the  contour  labeled  2  gives  the 

-2 

location  at  which  the  photon  flux  is  10   photons  per  unit 

area,  where  unit  area  is  one  extinction  length  squared. 

E.   FIELD  OF  VIEW  CONSIDERATIONS 

Figures  11  and  12  represent  the  relative  number  of 
photons  detected  at  a  receiver  located  zero  to  nine  and  36 
to  45  degrees  off  from  the  axis  of  propagation  at  one  and 
five  extinction  lengths  respectively  from  a  250  nm  laser 
source.  A  sixty  degree  (full  cone  angle)  field  of  view 
detector  would  receive  seventy-five  to  eighty  percent  of 
the  incident  radiation  for  the  0,  (zero  to  nine  degrees) 
and  the  0~  (36  to  45  degrees)  cases  at  one  extinction  length. 
At  five  extinction  lengths  these  percentages  tend  to  remain 
stable.   Thus,  given  a  minimum  required  photon  flux  at  a 
detector,  it  is  possible  to  select  the  proper  field 
of  view  for  that  detector.   Conversely,  utilizing  a  fixed 
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field  of  view  detector,  it  is  possible  to  determine  the 
maximum  effective  range  for  a  minimum  photon  flux.   (The 
results  for  280  and  300  nm  were  substantially  the  same 
as  the  250  nm  simulation) . 

F.   PULSE  SPREADING  CONSIDERATIONS 

Figures  13  and  14  represent  the  pulse  spreading  as 
received  by  a  detector  located  zero  to  nine  and  36  to  45 
degrees  off  from  the  axis  of  propagation  at  one  and  five 
extinction  lengths  from  a  250  nm  laser  source,  respectively 
The  pulse  spreading  function  is  a  function  of  the  position 
of  the  observer  and  time  dispersion  caused  by  multiple 
scattering.   The  relative  time  dispersion,  T(R,D)  is 
defined  as  follows: 


T(R,D)   =   |  -  1  (5) 


where  D  is  the  total  distance  traveled  by  each  photon  and 
R  is  the  distance  from  the  source  to  the  shell  penetrated 
by  each  photon. 

The  results  of  pulse  spreading  at  other  distances 
and  angles  was  examined  .    The   pulse   spreading 

function  is  nearly  constant  for  values  of  R  from  one  to  ten 
extinction  lengths,  thus  one  figure  suffices  as  an  approxi- 
mation to  any  distance  in  this  angular  range.   The  increase 
in  pulse  spreading  with  increasing  angle  is  illustrated 
with  one  additional  plot  at  36  to  4  5  degrees  off  from  the 
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axis  of  propagation.   Reproduction  of  the  results  at  each 
cone  angle  introduces  a  cumbersome  presentation  of  this 
information,  hence  only  one  representative  sample  is 
included. 

G .   SUMMARY 

The  absorption  and  scattering  coefficients  as  well  as 
the  phase  functions  from  previous  experiments  were  uti- 
lized and  the  behavior  of  multiple  scattered  radiation  was 
simulated  using  the  Monte  Carlo  method.   This  model  was 
used  to  study  the  characteristics  of  the  source  beam  at  a 
receiver  located  an  arbitrary  angle  off  from  the  axis  of 
propagation  at  a  given  distance  from  the  source.   With  such 
information  available,  the  geometric  pattern  of  possible 
communication  links  can  be  determined.   A  photon  density  as 
a  function  of  receiver  field  of  view  and  position  of  the 
observer  was  determined  as  well  as  the  effects  of  multiple 
scattering  on  pulse  spreading. 
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VI.   SIGNAL-TO-NOISE  CONSIDERATIONS 

A.  INTRODUCTION 

The  photon  flux  information  presented  in  the  previous 
section  may  be  utilized  to  estimate  the  signal-to-noise 
ratio  (S/N  or  SNR) .   This  information,  however,  must  be 
used  in  conjunction  with  the  characteristics  of  the  detec- 
tor.  Additional  considerations  therefore  include  the 
characteristics  of  the  optical  filter,  quantum  efficiency 
of  the  detector,  background  level  of  radiation,  etc.   This 
section  considers  the  characteristics  of  available  devices 
and  of  known  atmospheric  flux  levels  in  order  to  estimate 
the  flux  level  necessary  for  reliable  communication. 

B.  FUNDAMENTAL  RELATIONSHIPS 

The  detector  current  generated  by  a  light  signal  of 
power  Pq  and  of  frequency  v  is 

GneP 

La  =  -iet*  <6> 

where  G  is  the  gain  of  the  detector,  n  is  the  quantum 
efficiency  of  the  detector,  e  is  the  electronic  charge,  and 
h  is  Planck's  constant  [18]. 

The  dominant  noise  mechanism  for  detectors  operating 
in  the  visible  and  ultraviolet  region  of  the  spectrum  is 
shot  noise.  The  mean  square  of  the  noise  current  due  to 
this  noise  mechanism  is 
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<in2>  =   2eG2[ia  +  ^(PS  +  P  )]  B  (7) 


where  G  is  the  gain  of  the  detection  device  (typically  a 
photomultiplier) ,  i,  is  the  dark  current,  P   is  the  power 
of  the  background  radiation  present  due  to  sources  other 
than  that  from  the  desired  signal,  and  B  is  the  bandwidth 
of  the  receiver  [18,  25,  26,  27]. 

The  desired  expression  for  S/N  is  thus  the  ratio  of  the 
square  of  Equations  (6)  and  of  (7). 


2  nePq  2 

<i   >  (-t~) 

s/n  =  -V  »  % (8) 

<in  *  2e[id  +  Kv"  (ps  +  V]B 


In  several  specific  cases  the  dark  current  i,  is  negligible, 
whereupon  Equation  (8)  reduces  to  Equation  (9) . 


S/N  s  2S  'pT^pt'  (9) 


For  pulsed  signals  having  a  time  duration  of  x,  the  band- 
width B  is  effectively  1/2t.    The  relationship  for  S/N 
thus  reduces  to  the  number  of  photons  converted  to  electrons 
in  a  pulse,  degraded  by  the  ratio  of  the  signal  power  to 
total  power  received  over  the  duration  of  the  pulse. 

Each  of  the  expressions  for  S/N  may  be  modified  slightly 
by  consideration  of  a  duty  factor,  introduced  to  account 
for  the  fact  that  a  signal  is  not  constant  over  the  duration 
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of  a  pulse.   Consideration  of  this  introduces  a  multi- 
plicative factor  equal  to  Equation  (10) . 


Ds 

T/2 
/ 
.      "T/2 

m 

|        dt 

DN 

T/2 

/ 

-t/2 

VV°>J 

|        dt 

(10) 


where  D  /D   is  the  ratio  of  the  duty  factor  of  the  signal 

O    IN 

to  that  of  the  noise  over  the  duration  of  the  pulse,  x. 


C.   BACKGROUND  LEVELS  OF  RADIATION 

The  expressions  for  S/N  involve  critically  the  back- 
ground levels  of  radiation  due  to  other  sources.   In  communi- 
cation experiments  this  is  typically  due  to  scattered  sun- 
light.  Reliable  flux  levels  at  wavelengths  longer  than 
300  nm  are  reported  by  Valley  [28]  .   The  flux  levels  at 
shorter  wavelengths  are  a  result  of  both  direct  and  diffuse 
transmission.   These  levels  are  critically  dependent  on 
the  thickness  of  the  ozone  layer,  and  have  been  calculated 
for  various  atmospheric  conditions  [12] . 

The  photon  flux  calculated  by  Shettle  and  Green  is 
presented  in  Table  IV.   The  dramatic  drop  in  background 
radiation  levels  at  around  280  nm  is  evident.   The  signal 
that  would  be  transmitted  through  a  narrow  band  (10  nm) 

multilayer  dielectric  filter  at  300  nm  would  include  up  to 

13  2 

10    photons/cm   as  background  noise  photons,  whereas  at 
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wavelengths  below  280  nm  this  background  radiation  level 
may  be  considered  to  be  quite  small. 

TABLE  IV 

Photon  Flux  Due  to  Direct  and  Diffuse 
Transmittance  through  the  Earth's 
Atmosphere  [12] 

A  solar  angle  of  ninety  degrees  is  assumed,  along  with  a 

0.32  cm  ozone  thickness. 

2  2 

A  (nm)     Flux  (watts/m  nm)      Flux  (photons/cm  nm  sec) 


300 

1.07 

X 

10"2 

1.6 

X 

io12 

295 

5.15 

X 

lO"4 

7.8 

X 

io10 

290 

2.14 

X 

lO"6 

3.2 

X 

108 

285 

9.50 

X 

10"11 

1.4 

X 

104 

280 

9.09 

X 

io"19 

1.4 

X 

io'4 

D.   DETECTOR  CHARACTERISTICS 

The  detector  characteristics  which  relate  to  a  given 
S/N  expression  include  filter  characteristics,  quantum 
efficiency,  detector  area,  and  field  of  view.   Using  a 
combination  of  multilayer  dielectric  filters  and  absorption 
filters,  it  is  possible  to  construct  a  device  having  com- 
posite characteristics  which  will  pass  ten  percent  of  the 
incident  flux  at  26  5  nm  and  reject  the  ambient  solar  flux 
in  the  visible  to  a  degree  sufficient  to  reduce  the  overall 
background  count  level  to  approximately  ten  counts  per 
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second  [29].   The  transmission  characteristics  of  this 
filter  in  the  200  to  400  nm  region  are  indicated  in  Figure 
15.   It  is  thus  expected  that  somewhat  greater  than  one 
percent  but  less  than  ten  percent  of  the  photons  in  the 
250  to  280  nm  region  will  pass  through  the  filter.   Once 
the  photons  reach  the  photocathode  of  a  typical  photomulti- 
plier,  the  conversion  efficiency  is  usually  greater  than 
ten  percent  but  less  than  25  percent.   Thus,  for  purposes 
of  estimating  S/N,  the  combined  effects  of  filter  trans- 
mission characteristics  plus  quantum  efficiency  of  the 

photocathode  may  be  accounted  for  by  choosing  a  value  of 

-3  -2 

n  to  be  greater  than  10    but  less  than  10   .   With  unfor- 

seen  developments  it  is  possible  to  expect  a  combination 

having  a  composite  value  of  10   ,  however  such  a  device  is 

not  presently  available. 

The  area  of  the  detector  used  in  the  simulation  was 

chosen  to  be  nominally  one  square  inch,  typical  of  currently 

available  photomultiplier  tubes. 

E.   STATISTICAL  CONSIDERATIONS 

The  number  of  photons  converted  into  an  electrical 
signal  follows  Poisson  statistics: 


P(r)   =   e"n  nr/r!  (11) 


where  P(r)  is  the  probability  that  r  photons  are  observed 
in  a  sample  pulse  which  has  an  average  number  of  n  photons 
per  pulse. 
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FIGURE  15.   Transmission  Characteristics  for  a  Typical 
Composite  Filter 
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The  probability  of  observing  an  event  with  a  number  of 
photons  above  the  threshold  value  r,  .  is 

oo 

P(r>rth)   =     Z    e~n  nr/r!  (12) 

r=rth 

Graphs  of  this  function  are  readily  available  [30] .   Given 
an  average  number,  n,  equal  to  ten  and  a  threshold  level  of 
two  (to  distinguish  from  the  random  single  photon  events 
which  occur  at  the  detector  with  a  frequency  of  occurrence 
equal  to  ten  or  twenty  photons  per  second) ,  a  probability  of 
detection  of  such  pulses  is  .9995.   If  the  average  number 
n  is  five  with  the  threshold  number  equal  to  two,  the 
probability  of  detection  is  0.96.   Thus  an  average  value  of 
ten  photons  per  pulse  in  the  absence  of  background  radiation 
is  chosen  as  a  desired  level  to  avoid  excessive  error.   This 
means  that  a  desirable  value  of  S/N  is  ten. 

F.   SUMMARY  OF  SIGNAL-TO-NOISE  CONSIDERATIONS 

For  the  case  of  minimal  background  radiation,  Eq.  (13) 
gives  the  relationship  for  the  number  of  signal  photons 
required. 


V      1   S 

— ^—  =   -  (-)  (13) 

hv     n  n;  ^J; 


If  background  radiation  dominates  the  power  incident  on 
the  detector  (as  would  be  the  case  at  wavelengths  longer 
than  300  nm) ,  then  the  number  of  signal  photons  required  is 
given  by  Eq.  (14)  . 
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V        1   S    PBT   1/2 


Table  V  gives  a  summary  of  the  necessary  photon  flux 
external  to  the  detector  in  order  to  obtain  a  S/N  equal 
to  ten.   These  flux  levels  may  be  used  in  conjunction  with 
the  figures  of  Section  V  to  generate  the  geometric  patterns 
over  which  signal  transmission  may  be  expected. 

TABLE  V 

Photon  Flux  Required  in  order  to  Give 
a  S/N  Equal  to  Ten 

Assumed  parameters: 

Filter  characteristics  for  Figure  15 

Quantum  efficiency  of  ten  percent 

t   =   10   seconds  (significant  pulse  spreading) 

Area  of  detector  is  one  square  inch 

Note:    Filter  characteristics  at  300  nm  have  been  assumed 

> 

for  a  narrow  band  (ten  nm)  multilayer  dielectric 
filter. 


A  (nm) 

n 

PB 
hv 

.photons v 
sec 

V 

hv 

r 
-(photons) 

Required  flux 

2 

(photons/cm  ) 

250 

lO"3 



io4 

1.6  X  IO3 

265 

lO"2 



103 

1.6  X  IO2 

280 

lO"3 



104 

1.6  X  IO3 

300 

lO"2 

io13 

IO5 

1.6  X  IO4 
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VII.   DISCUSSION 

Although  not  currently  available,  three  milli joule 
laser  sources  of  different  wavelengths  (250,  280,  and  300 
nanometers)  capable  of  being  pulsed  at  an  appropriate  rate 
(200  pulses  per  second  for  2400  bit  per  second  data  rate 
using  12  bits  per  pulse  [4])  were  assumed.   (The  navy  has 
used  vocoders  with  2400  bits  per  second  capability  since  the 
1940's).   Table  V  lists  the  assumptions  leading  to  the 
required  photon  flux  in  order  to  achieve  a  S/N  equal  to 
ten.   By  utilizing  the  contours  of  equal  photon  flux  in 
Figures  8,  9,  and  10,  the  ranges  for  which  pulse  position 
modulation  communication  is  possible  were  predicted. 
Table  VI  lists  the  results  of  these  calculations. 

Predicted  performance  for  lasers  operating  at  250  and 
280  nm  are  comparable  for  sixty  degree  or  greater  off  axis 
detection.   For  detection  angles  less  than  sixty  degrees 
the  280  nm  laser  is  superior.   The  performance  of  the  300 
nm  laser  does  not  appear  to  be  degraded  by  the  background 
radiation  as  much  as  might  be  expected.   The  performance 
for  this  laser  from  zero  to  twenty  degrees  off  from  the 
axis  of  propagation  is  similar  to  the  280  nm  laser.   Beyond 
thirty  degrees  the  ranges  decrease  markedly,  indicating 
very  little  backscattering.   It  cannot  be  overstressed 
that  the  results  above  are  for  a  simulation  of  assumed 
laser  sources  with  assumed  atmospheric   propagation 
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characteristics  and  assumed  detector  characteristics. 
Clearly,  in  order  to  make  accurate  predictions  for  real 
systems  it  is  necessary  to  obtain  accurate  atmospheric 
characteristics  as  well  as  other  system  parameters. 
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VIII.   CONCLUSION 

Lasers  operating  in  the  middle  uv  do  not  appear  to  be 
suitable  for  long  range  communications;  however,  they  pre- 
sent definite  possibilities  for  short  range  (several  km) 
applications.   Since  radiation  in  this  region  is  absorbed 
exponentially,  rather  than  as  the  inverse  of  distance 
squared,  the  potential  for  covertness  is  highly  accent- 
uated.  Current  lasers,  however,  having  a  power  output  on 
the  order  of  0.1  joules  or  greater  with  low  pulse  rates 
(15  to  20  pulses  per  second) ,  are  not  suitable  for  voice 
communication  due  to  the  low  repetition  rates.   The  pulsed 
mode  of  operation  appears  to  hold  definite  advantages  over 
amplitude  modulation  of  a  continuous  wave  laser   from  simple 
signal-to-noise  considerations.   Good  possibilities  for 
successful  communication  are  offered  by  pulse  position 
modulation.   Millijoule  laser  sources  operating  with  200 
or  more  pulses  per  second  are  desirable  but  not  currently 
available. 

The  contours  of  equal  photon  flux  presented  in  this 
paper  are  useful  for  determining  the  communication  patterns 
of  possible  laser  sources.   The  pulse  spreading  and  field 
of  view  figures  are  likewise  useful  for  predicting  the 
characteristics  of  proposed  receivers.   A  multiple  scatter- 
ing model  is  definitely  required  for  proper  prediction  of 
photon  flux  over  expected  distances  of  communication. 
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APPENDIX  A 

MONTE  CARLO  SIMULATION  OF  THE 
MULTIPLE  SCATTERING  PROBLEM 

A.   FUNDAMENTAL  RELATIONSHIPS 

1.   Weighting  Factors  for  Monte  Carlo  Calculations 

The  following  calculations  served  as  a  model  for  the 
derivation  of  the  various  weighting  factors: 

PROBLEM:   Given  a  random  number  generator  that 
provides  numbers  within  a  specified  interval  along  the 
x-axis  (with  equal  probability  for  each  interval  dx) 
obtain  random  numbers  with  a  probability  p(y).   I.e., 
Find  y(x)  such  that  p(y)dy  =  dx. 

SOLUTION:   Let  Random  Variable  Y  be  distributed  in 
y   <_  y  <_   y2  and  its  density  function  be  given  by 

( p(y)  y0  <  y  <  y9 

p(y)  =  <        °      l     . 

'  Q         otherwise 
The  distribution  function,  P(y),  is  given  by 


y 

P(y)   =   P(Y  <  y)   =    /   p(y)  dy 


y 

Therefore,  P  (y)   =     /   p(y)  dy   ,  since  p(y)  =  0  for  y  <  y 

y~  c 


Now, 


p(y)  dy   =   dx 
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Integrating, 


Y(x)   =   yQ  +  P_1(Ax) 


where 


and 


y0  =  y(xQ) 

Ax   =   x  -  x 


This  concept  applied  to  specific  cases  yields  the 
following  fundamental  relationships  for  weighting  factors 
a.   Exponential 


P(y>   =  i  e  -^ 


(0  <  y  <  °°) 


y  (x)   =   -t  In  (1  -  Ax) 


b.   Henyey-Greenstein 


P(6)   = 


d-g2) 


2(l+g2-2g  cos  6)  3/2 


(0  <  0  <  tt) 


where 


g  =   <cos  0> 


9 (x)   =   cos 


-1 


iial- 

2g 


1 

2g 


_2g( 


i-g 
i+g  _ 
2g 


12 


Ax) 
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c.   Modified  Henyey-Greenstein  [Zachor  1977] 


p(6)   = 


_  i-g 


2  r 


i f  (3  cos  e-i: 

~    2    .       ox  3/2    2  ..   2.3/2 
(1+g  -2g  cos  6)         (1+g  ) 


0  <  9  <  7T 


where 


g   =   <cos  6>    and     f  is  a  weighting  factor. 


(x)   =   cos 


-1  )  1+g 
2g 


(g  -l) 


8g" 


[Ax-F(9)-  ^]2 


where 


F(9)   = 


£  d-g  ) 

4 


(1+gV77 


(cos  9  -  cos  9) 


d.   Uniform 


P(<1>)   = 


2tt 


0  <  4>  <  2tt 


4>(x)   =   2tt  Ax 


e.   Rayleigh 


P(9)   = 


3         2 
g-(l  +  cos  0) 


0  <  9  <  it 
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e (x)  =  cos  1(a  +  b) 


where 


B 


and 


b   =   8  Ax  -  4 

2 .   Position  and  Coordinates  in  Photon-Fixed  Coordinate 
System 

a.   Axis  Rotation 

Axis  rotation  is  accomplished  by  the  following 

transformation : 


x' 

y" 


cos  A6  cos  A(J)         -cos  AG  sin  Ac{>         -sin  A 9 

sin  A(f>  cos  A<£  0 

sin  A 9  cos  At}>        -sin  A9  sin  A<J)  cos  A9 


b.   Axis  Translation 

Axis  translation  is  accomplished  by  the  following 
transformation : 
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x" 

X' 

y" 

= 

Y' 

z" 

z  ' 

+   Ar 


3.   Angles  ( 9  ,  <j> ,  9  '  ,  <p   ')  for  a  Photon  Passing  Through 
a  Shell! 

As  As,  a* 

A  laser  is  located  at  the  origin  of  the  i  ,  j  ,  k 

o    o    o 

coordinate  system.   The  initial  direction  of  propagation  of 

the  collimated  beam  is  along  the  -k   axis.   6  is  the  angle 

between  the  position  vector  R  (pointing  from  the  location 

where  the  photon  penetrates  a  shell  toward  the  laser  source) 

and  the  vector  -k  .   (J)  is  the  angle  between  the  projection 

of  -R  onto  the  (i  ,j  )  plane  (denoted  by  -FL)  and  the 

o  Jo  P 

vector  i  . 
o 

A  photon  is  located  at  the  origin  of  the  V  ,  V  , 

f\  Ak  As,  As 

V   and  the  0,11,11   coordinate  systems.   As  the  photon 

z  x   y    z  2  ^ 

passes  through  the  shell,  01  specifies  the  angular  direction 


of  photon  travel  (along  U  )  relative  to  R  (the  line-of-sight 


position  vector) 


1  is  the  angle  between  the  projection 


of  U   (denoted  by  the  vector  W)  onto  the  (V  ,V  )  plane  and 

V  . 
x 

The  vectors  R. ,  R.,  and  R,  are  position  vectors 
i    j        k     r 

/\  A\  As 

in  the  U  ,  U  ,  U   coordinate  system  denoting  the  location 

x    y    z  ■*  3 

of  the  points  (1,0,0),  (0,1,0),  and  (0,0,1)  in  the  £  ,  j  , 


See  Figure  16. 
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^'M4 


FIGURE  16.   Diagram  of  the  Photon  Fixed  Coordinate 
System 
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k      coordinate   system.       (These   points    are   used    for   orientation 


purposes. ) 
a. 


=      cos 


'(- 


R    •    k 


o 


) 


cos 


-1 


R  k 


R   *     (Rk  -   R) 


b.         0   =      cos 


im 


<j>    =      cos 


-1 


p 


R"   •     (R.    -    R) 
l 


(o  <   e  <  tt) 


where 


V 


r2    -     [(RK    -    R)-R]2 


cos 


COS 


,  /u     .   r\ 

(hit) 


d.         4> '      =      cos 


,  /u     •   r\ 


(0    <    4>    <    2tt) 


(0    <    6'     <    tt) 


(0    <    <J>'    <    2tt) 


where 


/s  •v 


w      =      1   -    (V      •    U    ) 
z  z 


4 .   Mapping  Spherical  Coordinates  onto  Polar  Coordinates 
and  Coordinate  Transformation 

In  order  to  map  a  point  on  a  sphere  onto  a  circle (s) 
such  that  there  is  a  one-to-one  mapping  the  following 
transformations  were  utilized: 

Case  I.   0  <  9  <  tt/2 


r  =   R  /l  -  cos 


=   cos 


'[- 


Case  II.   tt/2  <  6  <  tt 


=   R  /l  +  cos 


cos 


-1 


-] 


In  order  to  plot  (0,<J>)  the  following  coordinate 
transformation  was  used: 

I.e.  -  Given  (9 ,  <J>)  ,  find  (x,y)  . 


R  /l  -  cos 


r   = 


R  /l  +  cos 


0  <_   0  <_  tt/2 
tt/2  <  0  <  tt 


x  =   r  cos  $ 


y   =   r  sin  <J> 
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5.   Phase  Functions 

A  factor  of  (yzr)  arises  from  the  normalization 
process  for  the   phase  function. 


2tt         it 
(i.e.         /  f      p(6)    sinGde    &$      =      1) 

0  0 


a.   Particulate 


PO)   = 


_   W 


4tt 


2  3/2 

(1+g  -2g  cos  9)  J/ 


f  3  cos  9-1) 

2  (iV)3/2 


b.   Rayleigh 


(1  +  cos2  9) 


167T 
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APPENDIX  B 

DEVELOPMENT  OF  THE  FUNDAMENTAL  RELATIONSHIP 
FOR  RAYLEIGH  SCATTERING 


p(6)   =   »>  (1  +  cos26)         0    9    7T 


e 

x  =    /   p(6)  d  cos 
0 


3        2 
;   -5-d  +  cos  6)  d  cos 

0    8 


3                           3  n      0 
,     Q    cos  6  x  i 
=  -  -g-  (cos  e  +  5 — )  i 


3 
3  ,    Q   cos   a    1 
x   =   -g-(cos  0  +  ^ — ^  +  2 


Now  solve  for  cos 


cos30  +  3  cos  0  +  (8x  -  4)   =   0 


Let 


y   =   cos  0 
b   =   8x  -  4 
a   =   3 
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y   +  ay  +  6   =   0 


(See  Reference  31  for  the  solution  to  a  cubic  equation.) 


or 


y   =   A  +  B 


or 


Now 


Y   =   " 


A  +  B     A  -  B 


/-3 


Y   =   " 


A  +  B     A  -  B 


/-3 


where 


and 


1     2    +1   4     27 


B   = 


f  2   T  4    27 


b2   a^ 

4  +  27 


2     3 
b     3 

4    27 


¥♦* 


V.2  V2 

Vi°;    V  +  x  >  ° 
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Therefore,  only  one  real  root  and  two  conjugate  imaginary 

roots  exist. 

Since  y  =   cos  0,  only  one  real  root  is  considered. 


6      =      cos"    (A  +   B) 


where 


and 


A  -  V-!  +  W+i 


'iiw?^ 


B      = 


b      =      8x   -    4 
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APPENDIX  C 
SAMPLE  CALCULATION  FOR  RANGE  PREDICTION 

PROBLEM:   Given  a  single  channel  280  nm  laser  utilizing 
pulse  position  with  power  output  per  pulse  equal  to  one 
milli joule,  determine  the  communication  pattern  for  a  S/N 
equal  to  ten.   Assume  the  detector  characteristics  presented 
in  Section  VI. 

SOLUTION:   From  Table  V,  the  minimum  photon  flux  for  a  S/N 

3  2 

of  ten  is  1.6  x  10   photons  per  cm  . 

Now, 


1  watt   =   5.0345  X  1018  A  Photons 

sec 


where 


A  =  .280  y  (A  is  measured  in  y) 
1  joule   =   1  watt  sec 


18 
=   1.40966  x  10    photons 


15 
1  millijoule   =   1.40966  x  10    photons 


-2 
From  Figure  9,  at  contour  2  the  relative  flux  is  10   of  the  initial 

photons  per  unit  area  where  the  unit  area  is  one  extinction 

length  squared.   From  Table  III  one  extinction  length  is 
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(1.0405  km   )    or  .9611  kilometers.   The  relative  flux  is 

-12    -2 

1.0826  x  10    cm     The  irradiance  is  given  by  the  product 

of  1.40966  x  10  5  and  1.0826  x  10~12  or  1526  photons  per 
square  cm.  By  measuring  the  distance  from  the  source  to 
the  detector  located  on  contour  2  (Figure  9)  and  recalling 

that  one  extinction  length  for  this  case  is  .9611  km,  the 

2 

ranges  for  an  irradiance  of  1526  photons/cm  may  be 

determined.   Table  VI  lists  the  results  of  this  calculation. 
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APPENDIX  D 
CHECKS  ON  POSSIBLE  ERRORS 

Each  computational  procedure  involved  in  writing  this 
program  was  verified  and  an  effort  was  made  to  ensure  the 
entire  program  functioned  properly.   The  procedures  uti- 
lized to  document  the  correct  behavior  of  this  routine  are 
described  in  this  appendix. 

The  various  weighting  functions  in  the  program  were 
verified  by  comparing  the  actual  distribution  of  photons 
generated  by  the  computer  to  the  analytical  solution  for 
the  distribution  of  random  numbers  obtained  using  a  program- 
mable calculator.   (The  statements  required  for  this  genera- 
tion remain  in  the  program,  but  the  output  of  the  results 
is  suppressed. ) 

Verification  of  axis  transformation  computations  was 
simple  and  straight- forward.   The  axes  were  rotated  about 
the  Z  axis,  the  Y  axis,  and  finally  about  both  axes  to 
ensure  proper  behavior.   (No  documentation  steps  of  this 
nature  have  been  left  in  the  program.)   The  output  of  9, 
<J>,0'/  and  <{> '  for  each  photon  at  each  shell  was  ordered  by  a 
computer  sort  routine  in  order  to  facilitate  the  preparation 
of  a  scatter  diagram  indicating  the  pattern  formed  by  the 
simulated  photons.   Inspection  of  such  a  pattern  did  in 
fact  reveal  the  presence  of  a  fault  which  was  corrected. 
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The  results  of  Duntley  [32]  in  describing  multiply 
scattered  radiation  in  the  propagation  of  laser  light  under 
water  are  most  useful.   Flux  as  a  function  of  cone  angle  is 
presented  for  distances  of  up  to  19  attenuation  lengths. 
The  scattering  phase  function  for  light  underwater  is  a 
sharply  peaked  function  which  may  be  approximated  with  a 
Henyey-Greenstein  function  with  a  g  value  equal  to  0.96 
[33].   Calculations  of  photon  flux  as  a  function  of  cone 

angle  were  carried  out  with  the  computer  program  developed. 

-3 
Total  attenuation  of  the  beam  down  to  a  level  of  10   could 

be  calculated  with  some  measure  of  accuracy  (only  10,000 
simulated  photons  were  utilized  for  these  calculations.) 
The  results  of  these  calculations  agreed  with  the  behavior 
observed  in  these  underwater  scattering  measurements,  simu- 
lating correctly  the  flux  levels  out  to  16  attenuation 
lengths  (the  limit  calculated) . 

Further  verification  of  the  accuracy  of  the  program 
was  obtained  by  comparison  with  calculations  by  Zachor  and 
Green  [24].   In  this  article,  photon  flux  as  a  function  of 
viewing  angle  was  calculated  for  several  distances  from  the 
source.   The  same  conditions  were  assumed  with  the  model 
used  here,  resulting  in  satisfactory  agreement  with  the 
information  published  by  Zachor  for  distances  of  1.43,  2.73, 
and  4.95  extinction  lengths  for  300  nm  radiation. 
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APPENDIX  E 
PROGRAM  DESCRIPTION  AND  DOCUMENTATION 

A.   COMPONENTS  OF  THE  PROGRAM 

1.   Weighting  Function  Subprograms 

RANEXP  is  a  function  subprogram  that  generates 
random  numbers  weighted  exponentially.  It  is  used  for 
distance  calculations. 

Function  subprogram  RANTH  generates  random  values 
for  theta  and  theta  prime.   Since  the  interaction  may  be 
either  Rayleigh  scattering  or  particulate  scattering, 
either  calculation  may  be  performed.   A  uniform  random 
number  (generated  by  RANDU)  within  the  interval  [0,1]  is 
obtained  and  compared  to  R   (the  ratio  of  particulate  to 
total  scattering) .   If  the  random  number  is  less  than  R 
the  particulate  scattering  calculation  is  made.   (Otherwise 
the  Rayleigh  calculation  is  performed.)   Due  to  the  com- 
plexity of  the  particulate  phase  function  it  is  difficult 
to  obtain  an  exact  solution;  therefore,  an  iterative  method 
is  employed.   An  exact  solution  is  available  when  f  is  set 
to  zero  and  a  more  efficient  calculation  is  utilized. 

RANPH  generates  a  random  value  of  phi  or  phi  prime 
uniformly  weighted  on  the  interval  [0,2tt].   The  steps  in 
this  function  subprogram  are  trivial. 
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2.   LITE  Subroutine 

The  main  subroutine  of  the  Monte  Carlo  simulation, 
called  LITE,  simulates  a  photon  which  randomly  interacts 
with  atmospheric  constituents.   (In  order  to  make  the  most 
efficient  use  of  computer  time,  non-interaction  was  not 
considered. )   Initially  the  photon  travels  along  the  -K. 
axis  (see  Figure  16,  Appendix  A),  a  distance  AR  (determined 
by  function  subprogram  RANEXP) .   A  uniform  random  variable 
in  the  interval  [0,1]  is  compared  to  R   (the  ratio  of  the 
total  scattering  coefficient,  Kq ,  to  the  total  extinction 
coefficient,  K) .   If  the  random  number  is  smaller  than 
R  ,  the  photon  is  scattered.   (Otherwise  it  is  absorbed  and 
the  photon  is  lost.)   If  it  is  scattered,  the  distance  of  the 
photon  from  the  laser  is  updated  and  a  new  value  of  R  is 
calculated.   Since  the  photon  has  been  scattered,  new  angles 
and  distances  must  be  computed.   This  is  accomplished  by 
generating  values  for  A9,  Acf),  and  AR  (using  RANTH,  RANPH, 
and  RANEXP,  respectively) .   The  photon  coordinate  system  is 
transformed  (translated  and  rotated)  to  a  new  coordinate 
system.   The  distance  from  the  source  is  calculated  and 
compared  with  the  radius  of  a  shell.   (A  series  of  concen- 
tric shells  with  radii  measured  in  extinction  lengths,  may 
be  formed  in  order  to  determine  the  photon  flux  or  angular 
distribution  of  photons  as  a  function  of  distance  from  the 
source.)   Various  calculations  are  then  made  to  determine 
whether  a  photon  has  penetrated  through  a  shell,  passed  in 
and  out  of  a  shell,  or  penetrated  an  outer  shell.   These 
calculations  are  necessary  to  ensure  proper  accounting  of 
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the  photon  flux.   When  a  photon  passes  through  a  shell,  the 
position  of  the  source  with  respect  to  the  penetration 
point  is  obtained.   Theta,  phi,  and  the  radius  of  the 
shell  (measured  in  extinction  lengths)  are  stored  in  a 
three  dimensional  array  for  future  use.   Theta  prime  and 
phi  prime  are  also  determined  as  a  function  of  distance 
(using  RANTH  and  RANPH) .   These  too,  are  stored  in  a  three 
dimensional  array  and  specify  the  direction  of  incident 
photon  flux  as  observed  by  a  receiver  located  a  given  dis- 
tance from  the  source  at  a  given  angle  off  from  the  axis  of 
propagation.   This  process  is  repeated  until  all  of  the 
photons  have  been  used  to  generate  information.   (It  should 
be  noted  that  many  extra  steps  incorporated  to  debug  the 
program,  remain  in  to  facilitate  future  implementation  by 
another  user. ) 

3.   DRLITE  Routine 

The  routine  DRLITE  is  used  to  drive  the  subroutine 
LITE.   Information  is  read  into  the  computer  and  the  results 
are  written  out  using  this  routine.   The  spatial  distribu- 
tion of  photons  as  a  function  of  position  and  detector  field 
of  view  is  calculated  as  well  as  the  total  number  of  photons 
in  each  shell.   The  output  of  other  available  information  is 
optional  and  is  controlled  by  IPRT  and  METHOD  statements. 
This  information  includes  pulse  spreading  (or  photon  time 
of  arrival)  as  a  function  of  position  and  field  of  view, 
the  relative  photon  flux  at  each  annular  ring  in  a  shell, 
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and  a  least  square  curve  fitting  of  a  modified  Henyey-Greenstein 
function  to  the  phase  function  information  generated  by  the 
computer. 

4 .   Miscellaneous  Subroutines  and  Function  Subprograms 
LSTHG,  LEAST,  GAUSS3,  and  RANDU  are  miscellaneous 
subroutines  used  in  the  program.   LSTHG  is  used  to  least 
squares  fit  a  modified  Henyey-Greenstein  function  to  the 
information  generated.   LEAST  and  GAUSS3  are  used  in  con- 
junction with  LSTHG.   EQN  is  a  function  subprogram  supporting 
LSTHG.   RANDU  is  a  uniform  random  number  generator  that 
provides  numbers  in  the  interval  [0,1].   It  is  located  in 
the  computer  library.   The  variable  IX  is  the  seed  and  may 
be  any  odd  integer  with  nine  digits  or  less.   It  may  not 
be  zero.   In  order  to  use  RANDU  the  following  FORTRAN 
statements  are  used: 

IX  =  948752759 

CALL  RANDU  (IX,IY,Y) 

IX  =  IY  . 

The  seed  IX  is  used  in  the  calculation  and  a  new  seed  IY  is 
generated  for  future  use.   A  uniform  random  number,  Y,  is 
in  the  interval  [0,1]  and  is  also  available  for  future 
calculations.   The  last  statement  (IX  =  IY)  resets  the  seed 
IX  to  a  new  value. 
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